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Abstract 

The behavior of an excess electron in a one, two and three dimensional classical liquid has 
been studied with the aid of Chandler, Singh and Richardson (CSR) theory [J. Chem. Phys. 81 
1975 (1984)] . The size or dispersion of the wavepacket associated with the solvated electron is 
very sensitive to the interaction between the electron and fluid atoms, and exhibits complicated 
behavior in its density dependence. The behavior is interpreted in terms of an interplay among 
four causes: the excluded volume effect due to solvent, the pair attractive interaction between the 
electron and a solvent atom, the thermal wavelength of the electron (A e ), a balance of the attractive 
interactions from different solvent atoms and the range of repulsive interaction between electron 
and solvent atom. Electron self-trapping behavior in all the dimensions has been studied for the 
same solvent-solvent and electron-solvent interaction potential and the results are presented for 
the same parameter in every dimension to show the comparison between the various dimensions. 
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I. INTRODUCTION 



The behavior of an excess eletron in a wide variety of fluids has been an interesting topic 
for many years0'Si'l@'illi'i0In the gas phase or the dilute liquid phase, the electron behaves 
almost like a free particle. As the solvent density increases, the electron exihibits different 
properties depending upon the nature of the solvent and electron-solvent interactions. At 
liquid density, the electron may become self-trapped in a cavity of solvent particles or remain 
quasifree depending upon the nature of electron-solvent interaction. Observed properties 
such as the electron mobilityS and the absorption spectra0 probe the nature of the electronic 
states in the fluid and phenomena of localization^. 

There are a broad range of problems in condensed matter physics that are intimately re- 
lated to the problem of excess electrons in deformable medium. These include charge trans- 
fer kinetics in biological reactions, metal-insulator transitions in fluids, polarons, phonon- 
assisted hopping of charge carriers in semiconductors and insulators, quantum-tunneling, 
etc. While the excess electron problem belongs to the general problem of electrons in dis- 
ordered materials, the liquid environment is in many ways different from the solid medium. 
In liquids, the constituent particles can diffuse, and local environment around the solute 
electron can be substantially different from that in solid. 

When an electron is solvated in a polar liquid such as water or ammonia, the strong 
anisotropic electron-solvent interaction causes significant local modification of the equilib- 
rium fluid structure&0 ? . The electron becomes localized in a small cavity because molecules 
in a solvation shell orient to create a potential minimum. Even simple fluids are found to 
exhibit electron mobilities that change by many orders of magnitude as the density of the 
fluid is altered slightly. In super critical helium, for example, the electron mobility drops by 
over 4 orders of magnitude as fluid density increased by a factor of 2 in low density regime! . 
The reason for this behavior is strong repulsion between electron and solvent atom. This 
causes the the depletion of the solvent atoms from the region of the electron and forms a 
highly localized state of the electron. 

In many other nonpolar fluids such as Ar, CH 4 , etc. the electron always remains in 
a state of high mobilitylll comparable to many semiconducting materials. An interesting 
density dependence of the mobility has been observed in them. It shows a minimum near 
critical fluid density and a maximum at liquid density. 
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Electronic states in reduced dimension are of considerable interest. For example, for a 
system less than two spatial dimension, electrons are localized with an infinitesimal amount 
of disorderlll. The interest in problems of electron or phonon propagation in a one dimen- 
sional random potential stems from the discovery and extensive experimental study of a 
certain class of organic or metallo-organic materials^. These matarials exhibit strongly 
anisotropic, quasi-one- dimensional behavior attributed to the fact that they consist of long 
chains, weakly interacting with each other. In many of these, the presense of a random po- 
tential has been proposed in order to explain their behavior. Electronic surface states play 
an important role in a wide variety of physical problems. For example, surface electrons on 
liquid helium has shown many interesting properties and led to important theoretical ad- 
vances such as the spectrum of bound electronic states, electron transport on the He surface, 
effects due to deformation of the He surface, and the possibility that the electrons may form 
a two-dimensional crystalEl in the field of low-dimensionality physics. 

The theory for the excess electrons in fluids developed by Chandler, Singh, and 
RichardsonS (CSR) is based on the path integral formulation of quantum theory which maps 
the behavior of the electron on to that of a classical isomorphic polymer^. The solvent- 
induced potential surface for the self-interaction of the isomorphic polymer is evaluated using 
an integral equation (e.g., reference interaction site model). With known potential surface, 
the polymer statistics is solved using variational approach^ that allows the determination of 
electronic properties and the structure of the liquid near the electron. The input of the the- 
ory is the pure solvent structure factor and the electron-solvent particle interaction potential. 
The CSR theory in its formulation is applicable to an adiabatic solvent (i.e., solvent particles 
are treated classically) , but has been extended to treat the effect of the quantum mechanical 
charge density fluctuation in the solvent particlesill. The calculated electron-absorption line 
shape and mobility are in good agreement with the simulations 7 and experiments!. The 
predictions of CSR theory were verified by computer simulations§Hl. 

For a one-dimensional system weil have shown recently that the repulsive and attractive 
parts of electron-solvent interaction potential lead separately to localization of electron, 
respectively, by creating cavity or forming a cluster of the fluids atoms around it. In two 
dimensional system weil have shown that dispersion of the wavepacket associated with the 
solvated electron is very sensitive to the interaction between the electron and the fluid atoms, 
and exhibits complicated behavior in its density dependence. CSR theory has been extended 
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to calculate the effective massEZI as well as the density matrix of the excess electron in fluidta. 

The CSR theory involves three or more characteristic lengths, depending upon the nature 
of electron-solvent interaction. These lengths are the thermal wavelength of excess electron 
A e = (^-)^ (where (3 is the inverse of temperature in the unit of the Boltzmann constant ks , 
m is the mass of a bare electron, and h is the Planck's constant divided by 2ir), characteristic 

length associated with the electron-solvent pseudopotential, and a length associated with the 

i 

mean volume occupied by each solvent atom, which is related to p* 17 , where p* = p s <7 D ,p s 
being the number density of the solvent and a is the diameter of the solvent atom, and D 
is the spatial dimensionality of the system. The behavior of the excess electron is expected 
to depend sensitively on these lengths. Laria and ChandleiS have attempted to explain 
the contrasting behaviors of the electron in super critical helium and xenon on the basis of 
different ranges of the electron-solvent repulsive interactions. 

In the present work we examine in detail the role played by different lengths and the 
spatial dimensionality of the system (in which we have considered the same solvent- solvent 
and electron-solvent model potential) to study the self-trapping behavior of the electron. 
The organization of the rest of the paper is as follows. In Sec. II we briefly review the CSR 
theory. In Sec. Ill we have presented the results and their discussions. Finally, Section 
IV presents concluding remarks. Appendix A provide some mathematical material for D- 
dimensional integration 



II. THEORY 



The system we consider is a single electron dissolved in a single component classical 
solvent. In the CSR theory, an excess electron is mapped, using a discretized version of the 
path integral formulation of quantum mechanics, onto a polymer of P interaction sites or 
beads. Under this isomorphism^, the electron can be viewed as a classical ring polymer. 

The total potential energy can be written as 

U=U ea (r,{R i }) + U a8 ({R i }) (1) 

with 

N 

U ss ({R i })= UssillU-Rjl), (2) 

i>j=i 
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and 

N 

Z7„(r,{R i }) = X;w«(|r-Ri|) (3) 

i=l 

Here, r denotes the position of the excess electron, Rj is the collection of the coordinates for 
a solvent atom, and N is the number of solvent atoms and u es (r) and u ss (r) are, respectively, 
electron-solvent atom and solvent atom-solvent atom interaction potentian. We consider a 
/^-dimensional system of spheres of diameter a in which the pair interaction between the 
solvent atoms is taken to be 

{oo for Irl < cr; 
n f I ^ 
for |r| > a. 

where |r| is the D-dimensional distance (for notational convenience the D dependence will 
not always be explicitly indicated). The electron-solvent atom interaction in a real system 
consists of a strong repulsion at short distance due to orthogonality requirements between 
wavefunctions of core electrons in the solvent particle and that of the excess electron and 
attraction at large distances due to dispersion interaction. However, in a system of neutral 
atoms the electronic states are determined primarily by the short range repulsive interac- 
tion or excluded- volume effect. The attractive interaction becomes important only at low 

densities. The interaction between the electron and solvent atom is taken to be 

{oo for r < d; 

- £ H£fcH> (oTr>d . < 5 > 
ar 

Here, d is the distance of closest approach between electron and solvent atom. 

In CSR theoryS the partition function Z for an electron in a bath of classical particles is 
written as the functional integral 

r r r 1 rP h 1 1 

Z = J Dv{u) J 4R,}exp —J du{-m|r(iO| 2 + MMR<})}-/^.(R<) ( 6 ) 

where r(u) is the electronic path in imaginary time which is periodic in time interval < 
u < j3h, i.e., r(0)=r(/3h) . To concentrate our attention on the electron degrees of freedom, 
the partition function given by Eq.(6) can be written as 

Z = Z S J Dr(u) ex V {-pS [r(u) - /3A^[r(«)]} (7) 

where Z s denotes the partition function of the solvent, A/i[r(u)] is the excess chemical 
potential for the fixed electronic path, 

1 rP % 1 

with, pS [r(u)} = - du-m\r(u)\ 2 (8) 
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and exp{— /3A/i[r(w)]} is called the influence functional which represents the solvent effects 
on the electon. In the continuum limiti, 

- A/i[r(«)] = Ps c es (0) + ^W~ 2 / du / du'v{\x{u) - r(u')l), (9) 

/ Jo Jo 

where c es (0) is the k = spatial Fourier transform of c es (r). 

v(\r(u) - r(u')\ = -J dr'J dv"c es {v', u) X ss{W ~ r"\)c es {\r" - r|,«') (10) 

Here r and u appear as independent coordinates, r is the distance between two sites and u 
measures the length along the contour of the polymer, and 

Xss(\r-r>\) = (5p s (r)5p s (r')} (11) 

is the density-density correlation function of the unpurturbed bath. In Eqs. (9) and (10), 
c es is the direct correlation function. Its value is determined from the equation^ 

p s h(r) = J dr'J dr"u(\r - r'\)c es (W - r"\) X ss(r") (12) 

where 

u(\t - r'|) = {phy 1 J d{u - u')u(\r - r'\;u - u'). (13) 

is the intrapolymer correlation function. 

Eq. (12) is solved for c es and h using suitable closure relational. Since all sites of a ring 
polymer on the average are equivalent, the site dependence disappears from Eq.(12) and 
only the zero- frequency component u;(|r|) of the equilibrium response function is required 
in Eq.(12). 

To complete the evaluation of excess chemical potential, the electronic path integral still 
has to be performed. Following FeynmamS and Chandler et. a/.,0 the excess chemical 
potential for the fixed electronic path is mimicked by a Gaussian functional, 

- pAp rei [r(u)] = -r o + -Wn)~ 2 / du / du'Y{u - v!) x \v{u) - r{u')\ 2 (14) 

2 jo Jo 

where T(u — u') is a solvent-induced force constant between different sites on the electron 
polymer and T merely determines the zero of energy. The Bogoliubov inequality provides 
a upper bound for the excess chemical potential, 

Ap < -/T 1 \nZ rc{ + < Ap[r(u)] - A/i rcf [r(u)] > rcf (15) 

6 



Here, Z Te { is the electronic partion function for the Gaussian reference system and < — > re f 
means the average over the reference system weight determined by S + A/i[r(w)]. Mini- 
mizing the right hand side of Eq.(14) provides the optimal Gaussian reference system. This 
procedure leads to the following equations. The correlation function for the intrapolymer 
correlation in k-space: 

u(k,r) = ex Ph^-] ( 16 ) 

where 

oo 

R 2 (r) = (\v(u) - r(u')\ 2 } = 4D £ A n [l - cos^r)] (17) 

n>l 

is the mean square displacement between two points on the electron path separated by a 
imaginary time increament < u — u' < (3fi with 

A n = ((3mQ 2 n + ln )- 1 (18) 

where Q n = and 

r /3h r d D kk 2 

7 „ = -{D(3hY l J o du[l - cos(O^)] J j^v(k) exp(-k 2 R 2 (u)/2D). (19) 

We solve Eq.(12) for h and c es using closure relation 

g = for r <d (20) 

c es = - ^('^ for r>d (21) 
ar 



We can express Eq.(12) and the closure (20) in the variational form 

(22) 
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where 

1 f d D k 

Irism = p s Ces(0) + -J j—^c 2 es (k)xss(k)Lu(k) (23) 

where c es (0) is the k = spatial Fourier transform of c es (r), D is the dimensionality of 
the space, dPk is the dimensionality dependence volume element, u labels the beads in 
the polymer ring (0 < u < (3K), m is the bare electron mass, and v(k) = —c 2 s {k)x ss {k) 



is the Fourier transform of the potential between beads, which is found in Eq.(lO). The 
information about the electron-solvent atom interation is contained in the closure relation. 

In Eq.(ll), u(r) is the intra polymer distribution function averaged over all beads of the 
ring polymer. In writing Eq.(ll) it has been assumed that for each polymer configuration, 
the solvent sees only average polymer rather than individual beads. The intra polymer dis- 
tribution function uj(r,r) is determined in the polaron approximatioiM,uj(k,u) is the Fourier 
transform of u(r,u) and is given by 

u(k,u) = exp[-k 2 R 2 (u)/2D]. (24) 



Eq.(22) is solved self-consistentlyEi for a given solvent and model potential representing the 
electron-fluid particle interaction. This solution gives information about v(r), w(r), R{u), 
7 n and g(r) [ = l+/i(r)]. Note that the quantity R(u) is the root mean square (RMS) value 
of the displacement between two points on the electron path separated by a time increment 
< u < /3h. The characteristic size or breadth of the polymer is measured by R(f3Ti/2). 
This is a measure of the spread of the wave packet associated with the particle. Since in the 
CSR theory a periodic boundary condition , r(0) = r(/3h), has been imposed on the path of 
the electron, R{u) is found to be symmetric about u = \(3Ti, i.e. it starts from zero at u = 0, 
attain a maximum value at u — \f3U and decreases for u > 4p reaching zero at u = (3%. g(r) 
gives information about the average packing of solvent particles around the electron. The 
variational parameter 7 n measures the strength of the electron fluid coupling. Quantities 
such as average kinetic energy, potential energy and effective mass etc., can be expressed in 
terms of 7 JEIzl as 



(K.E.) = ^k B T 



1+ ln 



(25) 



(3mVL 2 n + 7„ 

(P.E.) = p s Jdvu es (\r\)g(\r\). (26) 
-^ = 24£(4 7 rV + 7n A e 2 )- 1 . (27) 

" L n>0 

In this equation m* is the effective mass of the solvated electron. As mentioned earlier, 
we need two input for this theory. One is electron-solvent and another is the density- 
density correlation function of solvent which is related to the structure factor of the solvent. 

8 



For D- dimensional hard sphere solvent under cocsideration the Percus-Yevic (PY) equation 
can be solved analytically for D = llll and for D = 30. For D = 2 excellent results of 
thermodynamic and structural properties have been obtained by Baus and ColotH. In the 
above u, c es , and x, are the spatial Fourier transform of u, c es , and \- 

III. RESULTS AND DISCUSSIONS 

In presenting our results we mainly focus on the imaginary time correlation function 

R\hh) = (|r(~/3ft) - r(0)| 2 ) = AD £ A n (l - cosvrn) (28) 

1 1 n>l 

and the electron-solvent radial distribution function. Note that R{\(3K) gives a measure of 
the physical size of the electron chain or the spread of the wave packet associated with the 
electron. For a free particle, 

R{±0h) = A e . (29) 

In Fig.l we plot the reduced correlation length, S = R{\(3h) / \J D / A \ e which is the 
dispersion of the wavepacket associated with the solvated electron relative to the free particle 
in one, two and three dimension as a function of density for A e = 15a, a = a^ 1 and d/a 
= 0.29 for several values of the attractive interaction, f3e. From these figures we find that 
when the electron-solvent interaction is solely repulsive (/3e = 0.0), the electron is always 
gets trapped inside a solvent cage as solvent density is increased in one and two dimensions. 
In three dimension [see Fig. (1C)] the repulsive interaction is not strong enough (because 
d/a = 0.29) to localize the electron due to ordered structure formed when the solvent 
density is high. When the attractive interaction between electron and solvent atom is large 
((3e ~ 100), the reduced correlation length of the electron is very small at very low solvent 
density and it stays almost constant upto p* ~ 0.5. This behavior strongly indicates the 
electron is localized to a single atom irrespective of space dimensionality. As the space 
dimensionality increases, the range of the constant value of the solvent density decreases for 
example for D=l, S is constant upto p* ~ 0.6. One interesting feature we have noted from 
Figs. 1A , IB, and 1C is that as the attractive interaction increases, the reduced correlation 
length of the electron decreases as we increase the solvent density in the low-density regime, 
but the correlation turns upwardto have a maximum in the middle density regime, and 
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finally the electron is trapped in the high density regime. If we carefully look at Figs.l, we 
observe that as we increase the space dimensionality, we find these effect are pronounsed for 
large attractive interaction. The similar type of behavior has been observed for the other 
values of d/a. When d/a is less than 0.29 we found the reduced correlation length of the 
electron decreases as we increase the attractive interaction i.e. fie and is self trapped at 
lower density. Opposite trend is found for increasing d/a. 

From Figs. 1[A], 1[B] and 1[C] it is clear that at the begining when we start increasing the 
attractive interaction (fie ~ 20) the wave packet associated with the electron increases then 
further increase of electron solvent attractive interaction i.e. when fie > 50.0 the electron is 
localized on a single solvent atom. The reason behind the delocalized state at small vales of 
attractive interaction (fie) is the cancellation between the repulsive and attractive electron- 
solvent interaction. In Fig. 2 [A] to Fig.2[F] we plot the electron-solvent radial distribution 
function for one, two and three dimension at various solvent density and attractive iteraction 
(fie). Figs. 2 can be explained in consistent manner with the same physical picture given 
for Figs. [1] . When the attractive interaction is weak, the electron pushes the solvent atoms 
makes space to self-trap. When the attractive interaction is very large (fie > 100.0) the 
electron is trapped on a single atom irrespective of space dimensionality, as evidenced by a 
large peak in g(r) [see Figs. 2[F]. 

In Fig. 3 [A], 3[B], and 3[C] we plot the reduced correlation length as a function of X e /cr 
at fie = for D—l, 2, and 3 respectively. In one and two dimension ( Fig. 3 [A] and Fig.3[B] 
) we observed as we increase A e the electron is trapped in the solvent cage. If we compare 
Fig. 3A with Fig. 3B as the space dimensionality increases, the electron is strongly trapped 
in the solvent cage. One interesting feature we found in Fig. 3[C] at p* ~ 0.7 there is sharp 
transition from delocalized to self-trapped state. In less than two dimension when there 
is no attractive interaction, the electron will always be in a localized stateB but in three 
dimension there is always a tendendency the electron will be in delocalized state. On the 
other hand, as the temperature is decreased, the electron has tendency to be in the localized 
state. This competition between vaious length scales probably make the sharp transition 
around A e ~ 16a. 

In Fig. [4A], [4B] and [4C] we plot the reduced correlation length as a function of density 
(p*) for various values of d/a. We found for low d/a (for example d/a = 0.15) the attractive 
interaction dominates and the electron is self-trapped in the low density regime. On the 



10 



other hand, when d/a is large the self-trapping of the electron is dominated by repulsive 
interaction or self-trapping is by cage effect. In the case of electron in helium and xenon, 
the self-trapping of electron in helium and delocalized state of the electron in xenon has 
been explained on the basis of d/a value@ when elecron- solvent interaction is repulsive i.e. 
/3e = 0. 

IV. CONCLUDING REMARKS 

The CSR theory for the excess electrons in simple fluid ( consisting of spherical atoms ) 
has been studied in one, two, and three dimension. The detailed study led us to the following 
conclusions : 

1. The reduced correlation length, S, is very sensitive to the nature of the electron- 
fluid atom interaction, thermal wavelength of the excess electron ( A e ) and a length 

associated with the mean volume occupied by each fluid atom which is related to 

_ i 

p* D , where p* = pa D , p being the number density of the fluid atoms and a is the 
effective diameter of a fluid atom. The behavior is interpreted in terms of an interplay 
among the length scales noted above. When the electron-solvent attractive interaction 
is very large, the density dependence of the size of the electron polymer relative to 
the free particle is dominated essentially by the pair attractive interaction, and the 
electron is trapped in a single solvent atom irrespective of the dimensionality (D). On 
the contrary, if the attractive interaction is absent, the electron is trapped in a cage 
formed by the solvent atoms at higher density. When the attractive interaction is 
low (f3e ~ 10.0) due to cancellation between repulsive and attractive interaction, the 
electron is delocalized irrespective of the dimensionality. 

2. In the one dimensional case the reduced correlation length is almost independent of 
d/a when electron-solvent attractive interaction is absent. On the other hand, in two 
and three dimensions, S, the reduced correlation length of the electron is sensitive on 
the value of the d/a. 

3. In three dimension we found when d/a = 0.29 and /3e = 0.0, the temperature depen- 
dence of the reduced correlation length, S, shows a sharp transition (metal-insulator 
type) around X e /a ~ 15 at p* ~ 0.7. 
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4. The long range electron-solvent attraction dominates at low density, while hard core 
repulsion dominates at high densities. When both interactions are present, they can 
counterbalance each other. Consequently, at some intermediate density, the effective 
electron-solvent interaction can be quite small, resulting in the electron derealization. 

In the present work we have explored the electron self-trapping in simple liquid for various 
space dimension (D= 1,2, and 3) which have same model potential (solvent-solvent interac- 
tion potential and electron-solvent interation potential) in every space dimension. In every 
space dimension (i.e. D—l, 2, and 3) we have in the model system considered here there is 
repulsive as well as attractive electron solvent interation potential which is not considered 
in Refs. 25, 27, and 28. We have presented the dimensionality dependent result for self 
trapping behavior of electron in simple liquid. 
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APPENDIX A: D-SPACE DIMENSIONAL INTEGRATION 

In D-space we have for any function f(x) depending only on the distance x — |x|, 




(Al) 



or on x and one integration angle 0, 



poo p 

dxf(x,6) = S D ^ / dxx D ~ x \ d6sm D - 2 6f(x,9), 



Jo 



(A2) 



where Sd = DVd is the surface area of the unit sphere of volume 
_D/2),where T(l + z) — z\ is the T function. 
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